脉冲感应式电推进中平面型线圈激励等离子体的流动特征分析

doi:10.13675/j.cnki. tjjs. 180705

推进技术 ›› 2019, Vol. 40 ›› Issue (10): 2373-2382.DOI: 10.13675/j.cnki. tjjs. 180705

脉冲感应式电推进中平面型线圈激励等离子体的流动特征分析

成玉国¹,夏广庆^2,3

1.中国人民解放军91550部队41分队;2.大连理工大学工业装备结构分析国家重点实验室;3.大连理工大学辽宁省空天飞行器前沿技术重点实验室，辽宁大连;116024

发布日期:2021-08-15
作者简介:成玉国，博士，工程师，研究领域为空间推进技术。E-mail：hlcyg@126.com
基金资助:
国家自然科学基金 11675040 11702319国家自然科学基金（11675040；11702319）。

Analysis of Plasma Flow Characteristics under Excitation of Planar Coil in Pulsed Inductive Electric Propulsion

1.PLA 91550 Element 41,Dalian 116023,China;2.State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology,Dalian 116024,China;3.Key Laboratory of Advanced Technology for Aerospace Vehicles,Liaoning Province,Dalian University of Technology, Dalian 116024,China

Published:2021-08-15

摘要/Abstract

摘要： 感应等离子体推力器是一种具有较好应用前景的空间推进方式，受到复杂瞬态电磁场的影响，流场中热力学和动力学参数的时变特性难以实验测量。采用单流体磁流体力学理论，结合高温气体的热力学和输运模型，数值计算初始静止且无预电离的情况下，等离子体参数的二维分布，着重分析等离子体在前半周期内的流动特征。对内径r_i为0.05m，外径r_o为0.15m的线圈计算表明，其等离子体-线圈耦合距离约为0.032m，以Ar为工质情况下，有效冲量产生时，高价离子（Ar³⁺，Ar⁴⁺）与低价离子（Ar⁺，Ar²⁺）由于受电磁力影响不同，产生了局部分离电场，推动了电流片运动，且在电流片内部，前缘主要为低价离子，后缘主要为高价离子。单脉冲能量210.7J、峰值径向磁感应强度为0.5T、有效作用时间约12μs条件下，5μs时刻的Ar²⁺，Ar³⁺的最大数密度大于6×10²¹m^-3，且大于Ar⁺，Ar⁴⁺的数密度。电流片运动大于耦合距离后，受径向运动以及激励电流反转的影响，线圈表面等离子体磁场非线性特征显著，而前缘磁场维持规则分布。对外径为0.15m，0.3m和0.5m的线圈计算发现，r_o为0.5m时，电流片径向均匀长度达到0.3m，表明较大线圈尺寸除增加耦合距离外，可提高径向电流密度的均匀性。

关键词: 感应等离子体;高温热力学;电流片;流动特征;数值计算

Abstract: Pulsed inductive thruster is a promising plasma accelerator for space propulsion. The thermodynamic and kinetic parameters of the flow are difficult to measure in experiment due to complex transient electromagnetic field. In this work, the two-dimensional flow characteristics of the plasma, especially in the first half period, are numerically studied, by employing the single fluid MHD theory, in which the high temperature thermodynamic and transport models are included. In the simulation, the initial gas is static and without ionization. The calculations find that, for the excitation coil of r_i=0.05m and r_o=0.15m, the decoupling distance is about 0.032m and the local electric field driving the current sheet is generated by high level ionization ions (Ar³⁺,Ar⁴⁺) and low level ions (Ar⁺,Ar²⁺ ) caused by the difference in electromagnetic force when Ar is used as propellant. The front of the current sheet is comprised of low order ions, while the back is of high order ions. Maximum densities of Ar²⁺，Ar³⁺exceed 6×10²¹m^-3,which are greater than Ar⁺,Ar⁴⁺ at 5μs for single pulsed energy of 210.7J, peak radial intensity of 0.5T, period of about 12μs. As the current sheet moves far away from the coil, the magnetic field varies non-linearly on the coil surface and distributes regularly within the sheet. Analysis of the effect of coil dimension on the current uniformity suggests that, when r_o=0.5m, the radial uniform dimension can reach 0.3m, which is better than outer radius of 0.15m and 0.3m, which shows that as coil dimension increases, the uniformity of the current sheet in radial direction is enhanced, besides increasing the decoupling distance.

Key words: Inductive plasma;High temperature thermodynamics;Current sheet;Flow characteristics;Numerical calculation

成玉国,夏广庆. 脉冲感应式电推进中平面型线圈激励等离子体的流动特征分析[J]. 推进技术, 2019, 40(10): 2373-2382.

CHENG Yu-guo¹,XIA Guang-qing^2,3. Analysis of Plasma Flow Characteristics under Excitation of Planar Coil in Pulsed Inductive Electric Propulsion[J]. Journal of Propulsion Technology, 2019, 40(10): 2373-2382.

参考文献

[1] Mazouffre S, Grimaud L. Characteristics and Performances of a 100W Hall Thruster for Microspacecraft[J]. IEEE Transactions on Plasma Science， 2018， 46(2): 330-337.
[2] 李文博，胡俊锋，孙昊，等. 电推进空心阴极热子的寿命评估研究[J]. 推进技术， 2017， 38(9): 2146-2151.
[3] Martin A K. 2016 Performance Scaling of Inductive Pulsed Plasma Thrusters with Coil Angle and Pulse Rate[J]. Journal of Physics D:Applied Physics， 2016， 49(2).
[4] Polzin K A. Comprehensive Review of Planar Pulsed Inductive Plasma Thruster Research and Technology[J].Journal of Propulsion and Power， 2011， 27(3): 513-531.
[5] Polzin K A， Sankaran K， Ritchie A G， et al. Inductive Pulsed Plasma Thruster Model with Time-Evolution of Energy and State Properties[J]. Journal of Physics D:Applied Physics， 2013， 46(47).
[6] Lovberg R H， Dailey C L. Current Sheet Development in a Pulsed Inductive Accelerator[R]. AIAA89-2266.
[7] Dailey C L. Pulsed Electromagnetic Thruster[R]. AFRPL-TR-71-107.
[8] Russell D， Poylio J H， Goldstein W， et al. The Mark VI Pulsed Inductive Thruster[R]. AIAA2004-6054.
[9] Dailey C L， Lovberg R H. The PIT MKV Pulsed Inductive Thruster[R]. NASA-CR-191155.
[10] Mikellides P G. An Overview of the MHD Code， MACH[R]. AIAA2007-5614.
[11] Mikellides P G， Neilly C. Modeling and Performance Analysis of the Pulsed Inductive Thruster[J]. Journal of Propulsion and Power， 2007， 23(1): 51-58.
[12] Mikellides P G， Ratnayake N. Modeling of the Pulsed Inductive Thruster Operating with Ammonia Propellant[J]. Journal of Propulsion and Power， 2007， 23(4): 854-862.
[13] Mikellides P G， Villarreal J K. High Energy Pulsed Inductive Thruster Modeling Operating with Ammonia propellant[J]. Journal of Applied Physics， 2007， 102(10).
[14] Kim S， Soogab L， Kim K H. Wavenumber-Extended High-Order Oscillation Control Finite Volume Schemes for Multi-Dimensional Aeroacoustic Computations[J]. Journal of Computational Physics， 2008， 227(8): 4089-4122.
[15] Dedner A， Kemm F， Kroner D， et al. Hyperbolic Divergence Cleaning for the MHD Equations[J]. Journal of Computational Physics， 2002， 175(2): 645-673.
[16] Frese M H. A Two-Dimensional Magnetohydrodynamic Simulation Code for Complex Experimental Configurations[R]. AMRC-R-874.
[17] Murphy A B， Tam E. Thermodynamic Properties and Transport Coefficients of Arc Lamp Plasmas: Argon， Krypton and Xenon[J]. Journal of Physics D:Applied Physics， 2014， 47(29).
[18] Zhang X N， Li H P， Murphy A B， et al. A Numerical Model of Non-Equilibrium Thermal Plasmas. I. Transport Properties[J]. Physics of Plasmas， 2013， 20(3).
[19] Hirschfelder J O， Curtiss C F， Bird R B. Molecular Theory of Gases and Liquids[M]. New York: John Wiley， 1954.
[20] Colonna G， Laricchiuta A. General Number Algorithm for Classical Collision Integral Calculation[J]. Computer Physics Communications， 2008， 178(11): 809-816.
[21] O’Hara H， Smith F J. The Efficient Calculation of the Transport Properties of a Dilute Gas to a Prescribed Accuracy[J]. Journal of Computational Physics， 1970， 5(2): 328-344.
[22] Aziz R A， Slaman M J. The Repulsive Wall of the Ar-Ar Interatomic Potential Reexamined[J]. The Journal of Chemical Physics， 1990， 92(2):1030-1035.
[23] Rat V， André P， Aubreton J， et al. Transport Coefficients Including Diffusion in a Two-temperature Argon Plasma[J]. Journal of Physics D:Applied Physics， 2002， 35(10): 981-991.

脉冲感应式电推进中平面型线圈激励等离子体的流动特征分析

Analysis of Plasma Flow Characteristics under Excitation of Planar Coil in Pulsed Inductive Electric Propulsion

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 0

编辑推荐

Metrics

本文评价